WEAKLY NONLOCAL IRREVERSIBLE THERMODYNAMICS - 
THE GUYER-KRUMHANSL AND THE CAHN-HILLIARD 

EQUATIONS 



Abstract. Examples of irreversible thermodynamic theory of nonlocal phe- 
nomena are given, based on generalized entropy current. Thermodynamic 
currents and forces are identified to derive the Guyer-Krumhansl and Cahn- 
Hilliard equations. In the latter case Gurtin's rate dependent additional term 
is received through the thermodynamic approach. 



1. Introduction 

In the last decades there is a continuous interest in generahzed classical contin- 
uum theories to include nonlocal phenomena. The best treated and most popular 
type of nonlocal theories includes higher- and higher order space derivatives (gra- 
dients) of the relevant variables into the governing equations. Those investigations 
are frequently called as gradient or weakly nonlocal continuum theories. However, 
one cannot include higher order gradients in the governing equations at an arbi- 
trary place, all equations are to be compatible with basic physical principles, first 
of all with the Second Law of thermodynamics. More properly, we are to construct 
weakly nonlocal differential equations where the constitutive properties are deter- 
mined so that solutions of the governing dynamic equations result in nonnegative 
entropy production. 

Weakly nonlocal theories, related thermodynamics, are either supposing a gradi- 
ent dependent entropy function |^ , or are considering a generalized form of the 
entropy current § |, |. However, none of them uses the traditional methods of 
non-equilibrium thermodynamics, there are no thermodynamic forces and currents 
identified and no conductivity relations defined. In this paper we will find suitable 
thermodynamic forces and currents with the help of the generalized entropy current 
first used by Nyiri j^. Moreover, we will show that in case of heat conduction the 
classical Onsagerian method leads to the well-known weakly nonlocal generalization 
of the heat conduction equation called Guyer-Krumhansl equation. On the other 
hand, we can get an equation of the Cahn-Hilliard-type as an evolution equation 
of an extensive variable in the presence of an internal variable of nonlocal type. 
Several previous ad-hoc assumptions are explained by constitutive relations in the 
present treatment (e.g. the balance form equation of the heat current). 

1.1. Weakly nonlocal extended heat conduction. When treating heat con- 
duction phenomena in solid material, the starting point is the balance equation of 
internal energy 

(1) dtU + V-iq^O. 

Here u is the internal energy density and jg is the heat current density and 
dt denotes the partial time derivative. Extended thermodynamics considers the 
heat current density as an independent variable, introducing it into the basic state 



space and a non-equilibrium entropy is defined as a function on tliis extended state 
space. With, tfiese assumptions one can explain inertial, memory effects getting an 
appropriate phenomenological model of the wave heat conduction phenomena. The 
key to receive a well tractable model is to identify suitable thermodynamic currents 
and forces for a reasonable constitutive theory. This can be achieved by supposing 
that the non-equilibrium entropy function depends on the currents in a special way, 
as it was first suggested by Gyarmati 

(2) s(Mjg) = so(-u) - ijg • m • jg. 

Here m is the so called matrix of the thermodynamic inductivities. Further on 
we will assume that m is constant (as usual). The reciprocal absolute temperature 
is introduced as the derivative of the entropy function according to the internal 
energy H = ^ = 1/T. 

A natural physical assumption is that the entropy, being connected to internal 
energy, does not flow without the heat flow. Therefore, we accept, that the entropy 
current as a function of our given state variables (u, jg) has the property 0) = 0. 
Using Lagrange mean value theorem and considering this condition, one can get 
the next functional form 

(3) isiu,iq) ^B{u,iq) -iq, 

where the B called current intensity factor can be supposed to be continuous at a 
neighborhood of the local equilibrium state jq = 0. One can observe that B is one 
tensorial order higher than the corresponding current. After a short calculation, 
a consequence of our assumptions (|l|), (Q) and is, that the entropy production 
has the following form 

(4) a, = 9ts + V • = (B - l/TI) : V o j, + (V • B - m • dtjq) ■ j,. 

Here I is the second order identity tensor, : and o are the notations of the double 
contraction and the tensorial product, respectively. The entropy production is 
nonnegative and equals zero only in thermodynamic equilibrium according to the 
Second Law. 

We are looking for constitutive relations for B and dtiq that automatically fulfill 
the inequality above. The following force-current structure can be chosen 
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The linear approximation of Onsager can be written as follows: 
jq = Lii(V • B - m • + L12V o jq, 
(5) B-il = L2i(V-B-m-9tj,) + L22Voj,. 

Of course we could give more general relations among the thermodynamic cur- 
rents and forces, but this linear approximation seems to be quite general for our 
further purposes. From the equations above B can be easily eliminated and we can 
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give a general constitutive relation purely for 

jq = Lll(V^;^ - m • afjq) + LllV • ((L22 - L2lLj;i^Ll2)V O 

(6) + LsiLj^iij,) +Li2Voj,. 

This general relation gives the first approximation considering the entropy trig- 
gering effects of nonlocalities. The relations (^) and (|^) together result in the system 
of equations to be solved for u and jq. 

Let us now restrict ourselves to the isotropic case. According to the Curie prin- 
ciple (representation theorems of isotropic tensors), there is no cross coupling in 
the conductivity equations. In a strictly linear approximation (when L is constant) 
we get m = ml, Ln = 11, where m, I are a constant scalars and L22 can be written 
with indices as: 

{L22)ijki = hSikSji + hSiiSjk + hSijSki, 
Here li, I2, I3 and I are positive constant scalar material parameters to ensure 
nonnegative entropy production, m is also positive supposing concave entropy in 
the extended state space. Therefore the linear approximation of Onsager (^ reduces 
to 

jq = 1{V - B - mdtiq), 

(7) ^ ^ h{Vojg)+hiVojg)* -iql, 

where * denotes the transpose. Using this expression we can get 

(8) Imdtjq +jq = + lihAjq + [h + h)V O V • jg). 

This is exactly the Guyer-Krumhansl equation for the heat current, introduced to 
describe the thermal properties of some crystals at low temperatures Q . Originally 
it was derived with the use of kinetic physics but other derivations are based on two- 
fluid hydrodynamics . In the last decades there were several attempts to get these 
equations from pure non-equilibrium thermodynamics |l3|. However, 

all of these derivation contain several ad-hoc assumptions as supposing a balance- 
like dynamics for the heat current jg and specially simplified constitutive equations. 
Moreover, the derivations mentioned above does not have the heuristic power of 
irreversible thermodynamics, therefore the generalization of their assumptions in 
case of more difficult situations can be very ponderous if not impossible. 

In the special case of li = I2 — Iz = ^ we can get the Cattaneo-Vernote equation 
for the heat current. In this case the entropy current has its traditional form 
js — jq/T- Further, we get the Fourier heat conduction equation when m = 0, 
in the case of local equilibrium. Assuming the caloric state function in the form 
u = cT, the weakly nonlocal extended heat conduction equation can be given 
directly for the temperature as 

(9) ImcdttT + cdtT + lA^ + lc{h +I2 + h)AdtT = 0. 

Let us remark here that several generalizations are known for the classical en- 
tropy current expression from the investigations in extended thermodynamics |^ to 
the multifield theories that include finite length localization instabilities in damage 
mechanics models Q . The treated Nyfri form of the entropy current makes possible 
to exploit the Second Law in an Onsagerian spirit. 
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Another important remark is, that the heat current jg and the intensity factor 
B are internal variables in thermodynamic sense. If we do not eliminate B from 
the material equations, then it can be regarded as a current density of the heat 
current because the dynamic equation for the heat current derived above 

turns out to have a special balance form (more properly it will be an equation of 
Ginzburg-Landau type). Supposing that both the entropy and the entropy current 
depend on B, we can continue to introduce internal variables of nonlocal type in 
the following way: we construct a new current with the help of the given current 
intensity factor and introduce a new current intensity factor of B in the entropy 
current. In this way we can get a hierarchical approximation scheme of nonlocal 
phenomena based on thermodynamic considerations. 

2. Weakly nonlocal extensive variable - the Cahn-Hilliard equation 

Now, another example will be shown that similar structural assumptions of the 
entropy current can help to build up a current-force system in more general cases, 
too. Let us consider a continuum with a single scalar extensive variable a (e.g. 
mass density). Therefore a balance equation is given for a 

(10) 5ta + V-ja = 0. 

However, as we do not consider a simple local equilibrium approximation, we 
further extend the state space of extended thermodynamics and introduce an in- 
ternal variable ^ with the tensorial order of the current density. In this case our 
basic state space is spanned by the variables (a,ja,^)- We are to investigate the 
additional nonlocal terms in the governing equations, therefore, contrary to the 
previous example, we will neglect the memory effects, so the entropy function will 
depend only on the extensive quantities 

(11) s(aja,l) = so(a). 

The corresponding intensive quantity is denoted by T{a) — Let us inves- 
tigate now the functional form of the entropy current js(a,ja,^)- There are some 
natural physical assumptions that we apply: 

- There is no entropy flow in the absence of the flow of a and with zero internal 
variable: 

j,(a,0,0) =0. 

- In case of zero internal variable the entropy flow reduces to the classical form: 

Therefore, according to Lagrange's mean value theorem and the mentioned nat- 
ural physical assumptions we may write that 

js(a,ja,0 = rja + A-|. 

where A is a continuous constitutive function on the state space and its tensorial 
order is one order higher than the order of ^. We may observe that in this case 
the current intensity factor B connected to the heat current is replaced by its local 
equilibrium value F according to our assumption with neglecting the memory effects 
(expressed by the condition (|l^)). That simplifying assumptions on the entropy and 
the entropy flow yields the entropy production as 

(T, = r9ta + V • (Fja + A • I) = ja • Vr + A : V o ^ + (V • A) • ^. 
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Let us consider now the strictly linear approximation with constant coefficients 
in case of isotropic materials only. Now the choice of the currents and forces in 
straightforward. There are two coupled terms in the linear laws 

(12) = ZnVr + Zi2V-A, 

(13) ^ = Z2iVr + /22V- A, 

(14) A = hiW o ^) + h{V o ^)* + l^V ■ ^I. 

Here the matrix of the conductivity coefficients is positive definite to ensure 
nonnegative entropy production. We can eliminate ^ from the first two equations: 

^ = ^22^i2^ja + (^21 — ^22^i2^^1l)Vr- 

Eliminating A and ^ from the first equation, we get the following equality 

(15) ja =^llV(r-aAr)+Zi/22Aja + /22(/2 + Z3)VoV-j„, 

where a = {h + I2 + l3){hih2 — ^12^21) is positive. 



So we eliminated the internal variable and received (|lC|) and (15) as dynamic 
equations to be solved. 



Putting (|l5| ) into the balance (10) we finally obtain 



(16) dta + lnA{T - aAT) - ^22(^1 + ^2 + h)Adta = 0. 

This equation is astonishingly similar to the Cahn-Hilliard equation, but there 
are some essential differences. 

- The last extra term does not appear in the traditional equation. Gurtin 
received a similar additional term with the help of microforce balance (a 
strong additional assumption) in the special case when the variable a is the 
mass density [^5|. Supposing that ^22 = 0, this additional term vanishes. Let 
us remark that this rate dependent term appears in the linear Onsagerian 
conductivity equations at the same thermodynamic approximation level as 
all the others. 

- Under the second Laplace operator there is the intensive variable F, instead 
of the extensive a, contrary to the original Cahn-Hilliard equation. 

- All of the material parameters in (|l6|) are positive, as a consequence of the 
Second Law, a nonnegative entropy production. In the original equation the 
signs of the coefficients are fixed according to stability considerations that are 
not connected to the entropy balance. 

- The nonlinearities and anisotropics show a different structure than in the 
original equation. 

- In most of the previous derivations of the Cahn-Hilliard equations, a "gen- 
eralized chemical potential" appeared by analogy, and was inserted without 
any further ado into the mass balance. Here in our derivation we received 
a real generalization of the diffusion current in the same sense and with the 
same method as in the original irreversible thermodynamic approach. 

3. Discussion 

In this letter only some basic examples are treated. A more general theory 
of weakly nonlocal irreversible thermodynamics can be created according to the 
methods discussed above. It makes possible to incorporate and investigate several 
other classical nonlocal equations (phase field, Kardar-Parisi-Zhang, etc..) into 
the well established structure of non-equilibrium thermodynamics and investigate 
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them form the point of view of structural compatibiUty. The present treatment 
has several serious advantages if compared to other approaches. First of all the 
compatibility of those equations with the Second Law of thermodynamics can be 
investigated. These researches can deepen our understanding of the Second Law 
that seems to emerge as a general theory of material stability in a static and also 
in a dynamic sense . 

As we have already mentioned, the main advantage of the above scheme is the use 
of the force-current relations of Onsager, which offers a consequent approximation 
method. It is straightforward to go beyond the linear approximation and consider 
higher order terms in the series. Furthermore, with this scheme we can be sure 
that the amendment of the Second Law suggested by Muschik and Ehrentraut 
is fulfilled, that is the positivity of the entropy production is ensured by purely 
constitutive assumptions. 

One can say that the heuristic power of the Onsagerian approach can be a disad- 
vantage if contrasted with the more general exploitations of the entropy inequality. 
However, an application of the Liu procedure jl^, |9| shows that our assumptions 
on the entropy current are not only correct, but essentially one cannot construct 
more general constitutive relations that would satisfy the amendment of the Second 
Law. These problems and some other example applications are to be treated in a 
forthcoming paper |^ . 

Finally we remark, that there are not well established methods to derive weakly 
nonlocal equations in kinetic physics (see e.g. 0)- According to the present 
investigations, the role of the entropy current can be essential to develop a powerful 
approximation scheme, like the moment series expansion in case of memory effects 

i 
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